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Abstract 

We study four-dimensional simplicial gravity through numerical simulation with spe- 
cial attention to the existence of singular vertices, in the strong coupling phase, that 
are shared by abnormally large numbers of four-simplices. We attempt to cure this 
disease by adding a term to the action which suppresses such singular vertices. For a 
sufficiently large coefficient of the additional term, however, the phase transition dis- 
appears and the system is observed to be always in the branched polymer phase for 
any gravitational constant. 
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Nowadays, due to precise measurements of the gauge coupling constants at LEP, it is 
usual to consider the three interactions which are described quite well by the Standard Model 
to be unified into a Grand Unified Theory at about the energy scale of 10^^ GeV. On the 
other hand, it is well known that the gravitational interaction between particles becomes 
nonnegligible at the Planck scale, which is around 10^^ GeV. These two energy scales are 
remarkably close to each other considering the ambiguity involved in the derivation of the 
above values. The most natural interpretation of this fact is that, around such a high energy 
scale, all the four interactions, including gravity, are unified. There are two possibilities for 
such unified theories at present. One is string theory and the other is a unification within 
ordinary field theory including gravity. 

The problem we encounter when we try to formulate quantum gravity within ordinary 
field theory in four dimensions is that we cannot renormalize it perturbatively. In the 
path integral formalism, quantization of the geometry is performed by integrating over the 
metric field. If we use lattice regularization, which enables a nonperturbative study, general 
coordinate invariance is not manifest and whether it is restored in the continuum limit is a 
crucial problem. There are two kinds of lattice regularization of quantum gravity; namely 
dynamical triangulation and Regge calculus, of which the former is believed to restore general 
coordinate invariance in the continuum limit. It has been solved exactly in two dimensions 
[|l|] and its continuum limit is shown to reproduce Liouville theory 0, in which general 
coordinate invariance has been treated carefully. There is also a handwaving argument 
for the restoration of general coordinate invariance in the continuum limit of dynamical 
triangulation |^. 

Although four-dimensional dynamical triangulation seems to be difficult to solve analyt- 
ically, there is no potential barrier in studying it through numerical simulations. Employing 
the Einstein-Hilbert action as the lattice action and sweeping the gravitational constant, it 
has been discovered that the system undergoes a second order phase transition 0, which 
suggests the possibility of taking a continuum limit. The geometrical features of each phase 
have the following characteristics. In the weak coupling phase, the configurations are made of 
relatively small baby universes connected through pinches. This phase is called the branched 
polymer phase. In the strong coupling phase, on the other hand, the configuration has a 
mother universe with very high connectivity, which seems to be due to the existence of 
vertices shared by very large numbers of four-simplices. This phase is called the crumpled 
phase. 
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In the continuum theory, it is well-known that Euclidean Gravity may have a problem 
related to the unboundedness of the Einstein term 0. Decomposing the metric into the 
conformal mode and the transversal mode, one finds that the unboundedness comes from 
the fact that we can make the Einstein term arbitrarily negative by choosing rapidly varying 
conformal mode. This problem of the unboundedness of the Einstein term is, therefore, called 
'the conformal mode instability'. Since the coefficient of the Einstein term is proportional 
to the inverse of the gravitational constant, it is natural for us to relate this problem to the 
branched polymer structure observed in the weak coupling phase. 

In this respect, the strong coupling phase seems to be more promising than the weak 
coupling phase. At least, there is a mother universe, although it might be plagued with 
singular vertices shared by abnormally large numbers of four-simplices. Let us call the 
number of four-simplices sharing a vertex the vertex order. In this paper, we measure the 
vertex order distribution with high statistics. We observe two singular vertices which show 
themselves as a sharp peak. Moreover, we find that the peak moves towards a larger value of 
vertex order linearly as we increase the system size. Considering this fact to be a potential 
obstacle to taking a continuum limit, we add a term in the action which suppresses such 
singular vertices. For a sufficiently large coefficient of the additional term, however, we find 
that the phase transition disappears and that the system is always in the branched polymer 
phase. 

Let us first give a general outlook on this kind of approach. When we regularize four- 
dimensional quantum gravity with dynamical triangulation the integration over the metric 
is replaced with the random summation over all the four-dimensional simplicial manifolds. 
It is this point of the formalism that is essential in recovering general coordinate invariance 
in the continuum limit, as can be seen from the handwaving argument in Ref. 0. Note that 
the lattice action can, therefore, be chosen arbitrarily as long as it is local without being 
restricted by general coordinate invariance. We have to search for a sound second order 
phase transition, where we can take a sensible continuum limit, by changing the lattice 
action in various ways. Considering universality in quantum gravity as well as in ordinary 
field theories, we expect that the continuum theory, if it exists at all, does not depend on 
the way we construct it. 
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It is natural for us to start with the Einstein-Hilbert action 



S = I d'x^ (^-^^ 



as a first trial, where A is the cosmological constant and G is the gravitational constant. Let 
us denote the number of z-simplices in a simplicial manifold by Ni. One can easily find that, 
for a simplicial manifold. 



d^xy^ = cN, (2) 
J d^x^R = 2'kN2 - IQaNi, (3) 

where c is the volume of each four-simplex and a is the angle between two faces of a four- 
simplex, which is equal to arccos . Therefore the Einstein-Hilbert action can be 
expressed in terms of lattice variables as 

5lat = K4A^4 - «;2A^2, (4) 

where and K2 are related to A and G through, 

Ha = cA + -^ (5) 
27r 

^2 = ^. (6) 
The A^j's (z = 0, 1, ■ ■ ■ , 4) satisfy the following three relations. 

A^o - A^i + A^2 - A^3 + A^4 = X (7) 
2A^i - 3A^2 + 4A'3 - 5A'4 = (8) 
5N, = 2Ns, (9) 

where x is the Euler number, which is two for the spherical topology. Therefore only two 
of the Aj's are independent, which means that, actually, the lattice action is the most 
general one that can be written as a linear combination of Aj's. 

We consider an ensemble with a fixed A4 and with spherical topology. There are well 
established methods for generating such an ensemble through numerical simulations, and the 
technical details of our simulation shall be given elsewhere except for a few comments. 
The initial configuration is taken to be the surface of a five-simplex and we use the so-called 
(p,g)-moves ^ to update a configuration. Since the topology is not changed by the moves. 
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the spherical topology of the initial configuration is maintained throughout the run. Our 
code is written for arbitrary dimension following Ref. 

Let us turn to the results of our simulation. We first look at the second order phase tran- 
sition, which can be seen through thermodynamic quantities such as the average curvature 
per unit volume : 



(Rtot = J d'x,/^R) (10) 



_ Rtot 

No 

= 27r-/-10a. (11) 

Fig. |l| shows our results for (i?av) at various K2's. In contrast to the three-dimensional case 
0, no hysteresis has been observed. Also, one sees that the size dependence of the data 
changes abruptly at = 1-2 ~ 1.3. On the right there is little size dependence, whereas on 
the left, the curve goes lower and lower as we increase the system size. The derivative of the 
average curvature gives the susceptibility 

= 2"^;^ = — — ' ^^^^ 

which represents the fiuctuation of the total curvature. Fig. ^ shows the result. As is 
expected from Fig. the susceptibility has a peak around K2 = 1-2 ~ 1.3, which grows 
higher as the system size is increased. This implies that the correlation length of the local 
curvature diverges at the critical point |jTO|, where we may hope to take a continuum limit. 



These results are in reasonable agreement with the previous data in Refs. 0, |T^ . Since K2 
corresponds to the inverse of the gravitational constant, as is seen from (|^), we call the large 
K2 phase as the weak coupling phase and the small K2 phase as the strong coupling phase. 

However, as we mentioned earlier, the properties in both phases seem to be rather strange. 
In particular, considering the possibility of taking a continuum limit from the strong coupling 
phase, the singular vertices that are shared by abnormally large numbers of four-simplices 
might be an obstacle. We would like to study this phenomenon in detail. 

Let us look at the vertex order distribution, which can be defined as 

Pi^) = Tr(Y.^oiv),n), (13) 

where o{v) is the order of the vertex v. This quantity is measured every 100 sweeps and 
averaged over 100 configurations. In order to reduce the fiuctuations of the distribution, we 
smear the data over bins of size 10. 
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We first note that the average vertex order o{v) can be given as 

5A^4 



oiv) = ^Y.oiv) = —^. (14) 
In the last equahty, we used the relation 

Y.o{v) = 5N,, (15) 

V 

which comes from the fact that each four-simplex has five vertices. Thus one finds that the 



average curvature per unit volume -Rav can be written in terms of o{v) as 

No 

R,, = 27r-^-10a (16) 

iV4 

37r 



207ro('?;) - (10a - 4^) - — (17) 



~ 207vo{v) -0.614 (A^4>1), (18) 

where in the second equality we used the relation 

N2 = 2{No + N4-2). (19) 



Therefore, we can deduce from Fig. |l| how the average vertex order o{v) should behave. First 
of all, the average vertex order should be small in the weak coupling phase and large in the 
strong coupling phase. As for its behavior for increasing system size, in the weak coupling 
phase it should not change much, whereas in the strong couphng phase it should grow larger 
and larger. 

In Fig. 1^ we show the vertex order distribution p{n) for K2 = 0.0, 1.267 (near the critical 
point) and 2.0 with A^4 = 32000. For K2 = 2.0, the distribution damps quite rapidly for large 
vertex order and we have confirmed that it remains almost unchanged when we increase the 
system size. For K2 = 0.0, on the other hand, one finds that there is a peak of very large 
vertex order; as large as one third of the total four-simplices. We have also confirmed that 
the peak consists of two vertices. In Fig. ^ we show the size dependence of the vertex order 
distribution for K2 = 0.0. One finds that the very large vertex order grows linearly as one 
increases A^4. 

As was mentioned earlier, in the strong coupling phase, one already knows from the be- 
havior of (i?av) that the average vertex order should grow as the system size is increased. 
However, what we have found by looking at the vertex order distribution is that the dis- 
tribution does not shift to the right uniformly but that the two singular vertices move to 
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the right hnearly with increasing system size, leaving the continuum part which, as we also 
observe, moves to the right but much more slowly than the singular vertices. This seems to 
be a potential obstacle to taking a sensible continuum limit from the strong coupling phase. 

As we emphasized earlier, we still have the freedom to modify the action. We, therefore, 
attempt to cure the above disease by adding the following term to the action []. 

Su = ^T.io{v)-Cr. (20) 

We note that, since a change of the constant C can be absorbed into a redefinition of k,2 
and K4, we can fix it to an arbitrary value without loss of generality. We choose C = 5, 
considering that o{v) > 5 is guaranteed from the requirement that the configuration should 
be a simplicial manifold in the strictly mathematical sense. We have performed simulations 
with u = 5 ■ 10~^, 5 ■ 10~^ and 5 ■ 10~^. The vertex order distribution is measured every 100 
sweeps and averaged over 100 configurations as in the u = case. In order to reduce the 
fluctuations of the distribution, we smear the data over bins of size 5. 

In Fig. ^ we show the vertex order distribution for u = 5 ■ 10^^ with A^4 = 32000. One 
sees that there is no singular vertex with very large order for all K2. On the other hand, we 
find that for it = 5 • 10"^, the results resemble the ones for u = 0. 

Let us next turn to the curvature and the susceptibihty. Fig. |^ shows (-Rav) as a function 
of K2 for u = 5 ■ 10^^, 5 ■ 10~^ and 5 ■ 10~®. We have replotted the data for u = 0. For 
M = 5 ■ 10~^, the curves resemble the u = case. On the other hand, for u = 5 ■ 10^^ and 
5 ■ 10~^, the curves for A'4 = 8000 and 32000 coincide throughout the K2 region examined. 
Fig. 1^ shows the susceptibility xr as a function of K2 for u = 5 ■ 10~^, 5 ■ 10~^ and 5 ■ 10~^. 
Although we see a broad peak even ioi u> 5 ■ 10~^, it remains unchanged when we increase 
the system size in contrast to the u = case. 

These results imply that the second order phase transition which is observed in the u = 
case (Fig. ||) disappears for a sufficiently large u. We also observe that the configurations 
for u > 5 ■ 10~^ have much the same properties as the ones in the weak coupling phase of 
the u = case, i.e. no mother universe, low Hausdorff dimension, little size dependence in 
the vertex order distribution, etc.. Therefore we claim that in the presence of the additional 
term, with a sufficiently large coefficient in the action, the system is in the branched polymer 
phase for any K2 . 



For some studies on other types of modified action with different motivations, see Refs. 0,1 
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To summarize, we find that, in the strong couphng phase, there are two singular vertices 
with very large vertex order, which grows linearly as the system size is increased. We consider 
this to be a potential obstacle to taking a sensible continuum limit from the strong coupling 
phase. Adding a term like (|20|) turns out to have as strong an influence on the system as to 
erase the phase transition for a sufficiently large coefficient. A local term having such a large 
influence on the system is in itself an interesting fact. Also, that we find a branched polymer 
structure with a large average vertex order is worth mentioning. The fact that the size 
dependence completely disappears for a sufficiently large coefficient of the additional term 
provides an intuitive picture of the branched polymer structure : the typical configuration is 
composed of baby universes which are fluctuating independently, and increasing the system 
size merely contributes to increasing the number of baby universes without changing the 
local properties of the configuration. In order to take a sensible continuum limit, it seems 
that we need to fine-tune u or try other types of modification of the action 0]. 

We would like to thank H. Kawai, T. Yukawa and N. Tsuda for stimulating discussion. 
We are also grateful to B. Hanlon for carefully reading the manuscript. 
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Figure 1: The average curvature per unit volume (divided by a = arccosQj) is plotted 
against K2 for A^4 = 8000 (circles) and iV4 — 32000 (squares). 



Figure 2: The susceptibility (divided by a^) is plotted against K2 for — 8000 (circles) 
and Ni = 32000 (squares). 



Figure 3: The vertex order distribution for K2 — 0.0 (solid curve), 1.267 (dotted curve), 
and 2.0 (dashed curve) with = 32000. 



Figure 4: The vertex order distribution for K2 = 0.0 with A^4 = 8000 (dot-dashed curve), 
= 16000 (long dashed curve), = 32000 (dashed curve), N4 = 64000 (dotted curve), 
and A^4 = 128000 (sohd curve). 



Figure 5: The vertex order distribution for = 5 ■ 10 ^ with A^4 = 32000. The curves 
corresponds to K2 =-5.0, -3.0, -1.6, -1.2, -0.8, -0.4, 0.0, 0.8 and 1.6 respectively. 



Figure 6: The average curvature per unit volume (divided by a) is plotted against K2 for 
u = 5 ■ 10^^ (solid curve), 5 • 10~^ (dashed curve), 5 ■ 10^^ (thick curve) with A^4 =8000 and 
32000. The result for m = (thick dashed curve) with N4 =32000 is replotted. 



Figure 7: The susceptibility (divided by a^) is plotted against K2 ior u = 5 ■ 10 ^ (solid 
curve), 5 • 10"^ (dashed curve) and m = 5 • 10"^ (thick solid curve) with A^4 =8000 and 32000. 
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